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Classically chaotic systems relax to coarse grained states of equilibrium. This article reports work 
on quantized bounded relaxing systems, in particular the quasi-periodic fluctuations associated with 
the correlation between two density operators, whose classical limits as phase space partitions is 
evident. These fluctuations can distinguish classically chaotic and regular motions, thus providing a 
novel diagnostic devise of quantum chaos. We uncover several features of the relaxation fluctuations 
that are shared by disparate systems thus establishing restricted universality. The universal form of 
these fluctuations for generic quantized chaotic systems with time reversal symmetry is worked out 
within the framework of Random Matrix Theories and is found to be consistent with the numerical 
simulations performed on varied systems. 
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I. INTRODUCTION 



Classical dynamical systems have been classified into a hierarchy of deterministic randomness. We have 
well studied examples from integrable systems to Bernoulli systems, from the regular to those that are in a coarse 
grained manner identical to stochastic processes. An important notion in this context is that of mixing, as exemplified 
in the now classic coke-rum mixture of Arnold and Avez Q . Apart from the fundamental role it plays in deterministic 
randomness, alias chaos, it is the backbone of the foundations of classical statistical mechanics. While ergodicity is 
compatible with equilibrium, it is mixing that would ensure the drive to this state. The role of deterministic chaos 
in statistical mechanics is under active research, for instance it is not clear that thermal equilibrium is the result of 
only deterministic chaos or whether large system sizes are a neccessity. 

Low dimensional completely chaotic classical systems are however characterized by "equilibrium states" as the 
underlying invariant densities and mixing drives sufficiently well behaved initial states to this equilibrium. This 
affords the evaluation of averages much the same way as in equilibrium statistical mechanics. On the introdcution of 
quantum mechanics however, the picture is not quite so clear even for low dimensional systems. The phenomenon of 
quantum suppression of classical chaos in diffusive systems is now well studied, and several localization mechanisms 
have been put forward which inhibit quantum mixing pj , yet the issues in the situation of "hard chaos" for bounded 
chaotic systems relaxing to an equilibrium have not been sufficiently addressed. 

This paper studies quantum objects that are obviously connected to the twin issues of ergodicity and mixing in the 
classical limit. In the process we use models having a discrete spectra whose classical limit is known to be chaotic. 
The most convenient for our purposes is the quantization of two dimensional area preserving maps on the torus- a 
much studied subject. The Hilbert space is then finite dimensional and we have fully all the contradictions between 
quantum and classical chaos, while retaining the attractive feature of being able to do the quantum mechanics up to 
machine precision and size. 

Let the relevant classical phase space be ft. This could be for instance the energy shell on which the Hamiltonian 
flow is restricted. We may formulate mixing either in terms of phase space functions or densities. Let A and B be two 
subsets of the phase space such that they do not intersect. Let xa(x) and xb(x) be the two characteristic functions 
of these subregions. Let the flow integrated for a time t be denoted by /', this could be a continuous flow or a discrete 
map, we formulate this for the case of maps in which case t is an integer. The invariant measure of the flow is denoted 
by /it, in the present context of area preserving maps this is simply the physical area of the phase space region. 

The central quantity of interest is the correlation between the functions xa and xb- In more graphic language, as 
the subregion A evolves with the motion the correlation is the fractional area of its intersection with the region B. If 
in the long time limit this factorises into the fraction of the areas of A and B we have mixing. In other words the 
fraction of A systems in B is the fractional area of B in the long time limit. 

»(f(A)nB)/n(A)^»(B)/n(n) (i.i) 

Mixing systems are a step above ergodic ones in the hierarchy of classical dynamical systems. Ergodicity is the equality 
of the time average and phase space average of almost all points in the phase space. This can also be formulated as 
decorrelation on the average. Thus a system is ergodic if 

1 T 

lim - ^f\A) n B)/n(A) = n(B)/»(n). (1.2) 
t=l 

Classical mixing systems are characterized by a series of complex numbers or resonances that dictate the rate of 
decay of correlations. Thus for purely hyperbolic or Axiom A systems we can write 

fx(f(A) n B)/n{A) - n(B)/n(fl) =J2°i ex P (A 4 t). (1.3) 

i 

The Xi ( Ruelle resonances ) are in general complex numbers with negative real parts independent of the particular 
partitions A and B. Explicit calculation of these resonances is a challenging part of dynamical systems theory; there 
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are special models which are isomorphic to finite Markov processes for which these can be analytically found; for 
example the multi-baker maps . 

What we wish to study below are quantities in quantum mechanics that are naturally related to the classical 
correlation functions defined above. These are of course correlations between operators and when the dynamics is 
classically chaotic we expect the quantum correlations to also "decay" in some sense. Since we will consider systems 
with discrete spectra this does not strictly happen and eventually there are oscillations. We will characterize these 
oscillations and show that when the dynamics is classically chaotic they are well predicted by Random Matrix Theory 
(RMT). 

II. THE QUANTUM CORRELATION FUNCTION 

The relevant dynamical space on quantization is a Hilbert space ipv. For regions of phase space we consider 
subspaces of the Hilbert space. Since we have in mind classical maps on the two torus, let the Hilbert space be of 
finite dimensionality N. Let there be two distinct density operators on Hn (may or may not be projection operators), 
Pa and Pb- We will study the cases when these operators on Hn have obvious classical limits as functions on the 
phase space. We will restrict ourselves initially, and largely, to the case when the operators are projection operators 
such that Pa + Pb = In, where Jjv is the N dimensional identity matrix. Classically this corresponds to choosing 
the partitions A and B such that A U B = il. 

We sidestep the question of assigning an operator to a general subregion of phase space, and whether this is possible 
at all (say through the Wigner-Weyl transform of the characteristic function), by restricting this study to particularly 
simple subspaces (or density operators) both of and on Hm and going from the quantum to the classical instead 
of vice versa. Choosing an "arbitrary " subspace ( or operator) of Hjy obviously does not correspond to a proper 
subregion of the classical phase space and we will see that such subspaces do not provide interesting relaxation 
behaviors in that they do not distinguish between classically chaotic and regular motions. 

If we use pure states, say |n), we can construct the density operators as ^ |n)(n|. Let the dimensionality of the 
subspace A be tia = IaN; then tib = N — ha = JbN . The fractions Ja and fs determine the relative sizes of the 
partitions. We will mostly use projection operators that are diagonal in the position basis. In this case the classical 
region corresponding to 

Pa= J2 \ n )( n \ 

n=Q 

is the rectangle [0, Ja) x [0, 1) and the region corresponding to Pb is the rectangle [Ja, 1) x [0, 1), as the phase space 
is taken to be the unit torus (q,p) = [0, 1) x [0, 1). 

The dynamics is an unitary operator U acting on the states of Pat. The central quantity of interest is then 

C AB {t) = Tr{U t P A U- t P B )/N, (2.2) 

which gives us the overlap of the subspace A propagated for time t with the subspace B. We may note that this is an 
important and natural physical quantity to study and has appeared before in several contexts, for example ref. It 
can also be viewed as an analogue of Landauer conductance for bound systems and as such should be of considerable 
importance in quantum transport. If in fact there is factorization in time then Cab(oo) must be compared with 
Tr(P A )Tr(P B )/iV 2 = f A f B . 

Since U has a discrete and finite spectrum, the correlation can only be a finite sum of purely oscillatory terms 
and can therefore display decay only over short time scales (of the order of Heisenberg time). The decorrelation is 
therefore not expected. As an extreme example if we consider the case when the projection operators are constructed 
out of the basis functions of U the correlation is zero for all time. This is an extreme case and for instance it is 
not clear what the classical limit is of the Wigner-Weyl transform of projection operators constructed out of such a 
basis. We will consider "generic" subspaces and projection operators and any claim of universality made below has to 
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be viewed with this caveat. We may note that this plagues Random matrix theory descriptions of eigenfunctions as 
well, as it involves basis dependent quantities. Here the basis dependence enters as the basis in which the projection 
operators are diagonal. 

It was suggested earlier using quantized multi-baker maps that the quantum correlation function approached the 
classical correlation function in the classical limit corresponding here to N — ► oo 0. Certain remarkable features of 
quantum relaxation were noted there including relaxation localization and effects of symmetries on transport. Here 
we wish to study the fluctuation properties more closely and uncover possible universal features. Therefore we study 
rather standard models like the Taylor-Chirikov (standard) map the kicked Harper system ]l(J and the baker 
map |ll[ qalthough most of these models do not have explicit expressions for the classical Ruelle resonances and arc 
probably not Axiom A systems (except the baker map) . However it is very reasonable to assume that in certain well 
known parameter regimes they can be for all practical purposes mixing systems. 

Normalizing the correlation so that the classical limit, if it exists, is for large times unity we study the quantity 
c(t) = Cab (t) / Ja/bj where from now we will acknowledge implicitly the dependence of the correlation on the 
particular choice of partitions. Fig. 1 shows an example of the generic behaviour of c(t), this particular data being 
for the standard map at two different values of the inverse Planck constant N for the same classical value of the chaos 
parameter. The initial relaxing behaviour is not clear from the figures as the time scales are much larger than the 
inverse of the principal Ruelle resonance. The projection operators used are diagonal in the discrete position basis, Pa 
is given by Eq.( 2.1) and Pb — In — Pa- The position eigenbasis is denoted as |n) and in this figure $a = /g = 1/2. 

We have used dimensionless scaled position and momentum co-ordinates and time is measured as multiples of the 
period of the kick, taken as unity. Modulo one conditions restrict the phase space to the unit torus so that the 
dimension of the Hilbert space N is related to the Planck constant as N = 1/h and the classical limit corresponds to 
N — > oo while the parameters of the map such as the kick strengths are scaled and dimensionless. 

The first observation is that the quantum correlation is in fact quite close to the classical value of unity, second 
is that the average of the oscillations in relation to unity will give us some information about whether quantum 
mechanics is inhibiting transport or otherwise, thirdly this average must in the classical limit tend to unity, fourth is 
the observation that the fluctuations are getting smaller in the classical limit and must tend to zero. Some of these 
observations have been made and substantiated earlier (TjJ^], here we will elaborate and present more results. 

Previous work on the use of time developing states to study the quantum manifestations of chaos have basically 
concentrated on the survival probability of a pure state ||. This corresponds to the choice Pa = Pb = IV'KV'I with 
an arbitrary normalized initial state \ip) and thus j a = Jb = 1/N. Such survival probabilities averaged over random 
initial states show a distinct difference between chaotic and regular systems. In contrast what we study in this paper 
are density (or projection) operators and not pure states. It becomes quite important that /aN and fsN are large 
and are of the order of the Hilbert space dimensionality N itself. Besides any "arbitrary" state can be chosen for 
the calculation of the survival probability while we will restrict ourselves to those projection operators that can be 
interpreted as phase space regions in the classical limit. This is to ensure that the classical transition to chaos is fully 
reflected in the quantum relaxation, as will be illustrated below. 



III. THE FLUCTUATION PROPERTIES 

The quantum "equilibrium" as opposed to the classical is an highly oscillatory state, with the fluctuations coming 
from the discrete nature of the spectrum. We will denote the time average of the correlation function by < c > and 
its variance by a 2 . These quantities can be written in terms of the eigenfunctions of the evolution operator U, which 
satisfy the following equation. 

U\ij> m ) = e lE ™\4> m ) m = 0,...N-l (3.1) 

The real numbers E m are the eigenangles of the quantum map, and we are assuming that there are no exact degenere- 
cies as is the generic case. The correlation is then 
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c(t) 



1 



Nf A (l - f A ) 



Nf A -l N-l 

E E E e^-^^^I^X^J^Xml^)^!^) 



(3.2) 



mi,m 2 rii=0 n 2 =NfA 

Thus we can express the average quite simply as 



< c>= 



1 



N—I/aN—1 N-l 



E E E IW^HfaU 



N/a(1 Ja) m=Q n _ n „._ /vA 



(3.3) 



Thus the average is a measure of the distribution of the eigenvectors in the subspaces A and B. We suggest that 
< c > < 1 generically, indicating a certain reluctance to participate equally in both the partitions in proportion to 
their sizes. A quantity of natural interest is the variance of the fluctuations which exist irrespective of the symmetries 
of the system. This can also be expressed in terms of the eigenfunctions of the system. In the particular case when 
Pa + Pb — In, writing for j a simply /, we get after some simplifications 



a 2 = 



N 2p {l _ f) 2 



E 

m 2 >mi 



fN-1 
E W^K^mJ 

71=0 



(3.4) 



Thus the variance measures a correlation between distinct pairs of eigenfunctions. There is a coherent partial sum 
over Nf states that makes the variance non-trivial. 



A. The Models 

Apart from time reversal symmetry there may be phase space symmetries that are either preserved or not in the 
quantum mechanical models. The results in ||] are for the case when there is an additional phase space symmtery 
present which is reflected quantally as a parity symmtery of the eigenfunctions. This being a special (but important) 
case we continue our study here for the more generic case when the eigenfunctions possess no symmetries. Then the 
only important symmtery is that of time reversal. 

The Taylor-Chrikov, or the standard map, or the kicked rotor, used in this paper is described below |§,|l5|. Let the 
classical map be (qt+i = It + Pt+i, Pt+i — Pt — V(qt+i)), both q and p taken mod 1. V' is the derivative of the 
kicking potential which is assumed to have unit periodicity. The toral states are assumed to satisfy certain boundary 
conditions specified by a point on the dual torus. Let \q n ) and \p m ) be the position and momentum states then 
\Phi+n) — e~ 2ma \p m ) and |g n +jv) — e 27Tlb \q n ), where (a, b) are real numbers between and 1; N is the dimensionality 
of the Hilbert space. If b = 0, upon canonical quantization we get the finite unitary operator 

E/ n , n , = B -y= eM-^iNVi-^j) exp(*-(n - n'f). (3.5) 



We have used for the standard map V(q) — K cos{2nq) / {2n), and a — 1/2 makes the quantum map possess an exact 
parity symmetry about q = 1/2. Thus we break the symmetry by simply choosing different boundary conditions for 
the wavef unctions. If b =/= then time reversal symmtery is also broken. Therefore there are two classes of models 
with no special wavefunction symmetries. We call the case (a = 1/2, b — 0) the symmteric standard map. The values 
of N are restricted to the even integers. 

The kicked Harper map ]l0| on the torus is similar to the above system except for the momentum dependence. If 
the Hamiltonian is 

oo 

H = V 1 (p) + V 2 (q) S(t-n), 

n— — co 

the quantum map is 
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Un, n > = ^exp(-27riJVV a (^)) £ exp(-2«iVy 1 (^))exp(^(m + 6)(n-n')) (3.6) 

m=0 

For the map used below we have taken Vi (p) = —31 cos(27rp)/(27r) and V^(g) = — 32 cos(27r<7)/(27r). If (a = 6 = 1/2), 
this once again ensures special symmetry properties of the quantum map, a case we call below the symmetric Harper 
map. We can retain time reversal and break parity by choosing 6=1/2 and a ^= 1/2, and or break both symmetries. 
The classical map is (qt+i = qt +9i s in(27rp t+1 ), p t+1 = p t — g 2 sin(27rg t )), again the mod 1 rule is assumed. If g\^ are 
equal, as assumed in this paper, the transition to chaos occurs around 31.2 = 0.63. In both the kicked Harper model 
and the rotor the time between kicks has been taken as unity, as anyway there are two parameters in the quantum 
problem, the scaled Planck constant N and the kick strength (K or g). 

The bakers map fll|]l2|| comes in several varieties. It has been shown that the usual bakers map ((qt+i = 2qt, pt+i = 
Pt/2) if < q t < 1/2 and (q t+1 = 2q t - 1, p t +i = {pt + l)/2) if 1/2 < q t < 1) is isomorphic to the (1/2, 1/2) 
Bernoulli process. The quantum bakers map on the torus is then the unitary operator 

u ' G »( a 7 ol,)- <"» 

where Gjv is the finite Fourier transform matrix with elements 

(G N ) m , n = -±= exp(-27ri(m + l/2)(n + l/2)/N). 
v N 



We have assumed anti-periodic boundary conditions (a = b = 1/2) and this is known to preserve classical symmetries 
of the usual bakers map. 

The generalized bakers map used in this paper is a dynamical system implementing the (2/3, 1/3) Bernoulli process 
and is the classical map ((qt+i = 3q t /2, p t +i = 2p t /3) if < q t < 2/3 and (qt+\ = 3q t - 2, p t +i = (pt + 
2)/3) if 2/3 < qt < 1). The quantum map is the unitary operator 

U = G N 1 ( G ™' 3 ^ ) . (3.8) 



G 



N/3 



IV. NUMERICAL RESULTS 

Time dependent quantum chaotic systems such as the kicked rotor on the cylinder are known to suppress classical 
chaos and lower diffusion Q . The corresponding suppression in the case of bounded systems could be the average 
relaxation such as measured by < c >. We note that the sum over n and n' expressing the average can be factored into 
two single sums, and we see that if for instance we had a parity symmetry forcing the wavefunction to be essentially 
identical in the subspaces A and B we would have < c >= 1, which is indeed the case for the fluctuations shown in 
Fig. 1. The average < c > could also be greater than unity in the presence of symmetries, or if the partition B is 
identical to A. 

Fig. 2a shows the average for the quantized baker map implementing the Bernoulli scheme (2/3, 1/3) as a function 
of the inverse Planck constant N . The partitions are such that /a = /s = 1/2. Shown is the deviation of the average 
from unity and an approximate power law is observed for large N. Thus we can write 

< c>~ 1 - aN- 1 , (4.1) 

and in this case 7 w 3/4. The relaxation localization is implied by a being positive. Fig. 2b shows the same quantity 
for the symmetric standard map and we find for the partition j a — 1/4, Jb = 3/4 that 7 ~ 1, with a different value of 
a. We note that in both the cases there is a small oscillation about the fitted straight line. Further results, especially 
on the variation of the average with the parameters of chaos, are in ref. B. 
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Fig. 3 shows the standard deviation er scaled to aN as a function of the kick strength in the case of the symmetric 
standard map. The transition to classical chaos at K rs 1 is visible in the variance of the fluctuations as a point 
at which it attains an approximately constant value. The relaxation fluctuations are therefore significantly different 
depending on the dynamical nature of the classical limit, being in general larger for regular systems than chaotic ones, 
and thus provide a novel diagnostic devise for the study of quantized chaotic systems. Further results may be found 
in ref. |^| . Here we will turn to a more closer appraisal of the ingredients of the variance itself. 

In a slight modification of the equation for the variance, define 



S = 2 £ 



m 2 >mi 



no+l— 1 



^2 M^mj.XV'mil 
ti—uq 



(4.2) 



The quantity S is simply related to the variance of the fluctuations. The relevant partitions have been generalized 



slightly so that Pa — Xm^=n _1 \ n )( n \ anc ^ classical partition is the vertical rectangle of length l/N whose origin is 
at rio/N. The B partition is the complementary space. S is potentially a function oil, N and no. For generic systems 
the origin no should not play a role and we will expect this dependence to drop out. In fact the case of the symmetric 
standard map when the partition is constructed from the position basis is a counter-example to this statement. But 
what we will go on to establish is that for generic quantum chaotic systems and generic partitions, S is asymptotically 
(large N) only a function of / = l/N and that the functional dependence is determined by time reversal symmetry. 

In Fig. 4a is plotted the behaviour of S for the symmetric standard map and the dependence on no is clearly 
visisble as the non-analyticity as a function of the fraction /. On the other hand Fig. 4b illustrates a smooth 
functional dependence independent of no when the parity symmetry is broken by choosing a = 0.35 instead of 0.5. 
The universality of this curve has been checked against other models such as the kicked Harper map. In the next 
section a form of this smooth function is derived within the framework of RMT. Fig. 5 consolidates three cases of the 
kicked Harper map. The function S with both parity (R- symmetry) and time-reversal, with just time-reversal and 
the case when there is neither parity nor time-reversal symmetry. The fluctuations clearly distinguish these situations 
and illustrates the intuitive expectation that the fewer the symmetries the smaller would be the fluctuations about 
the "quantum equilibrium" . 



V. FLUCTUATIONS ACCORDING TO RMT 

Random Matrix Theory []l3|| has been quite successful in the description of many aspects of complex quantum 
systems. It is natural to expect that the fluctuations be also described adequately within this theory. Already work 
on the survival probability has been dealt using RMT with success. The fluctuations studied above are characterized 
by certain higher order correlations between eigenfunctions and the evaluation of these within the RMT framework 
for the case of time reversal symmetric generic systems is dealt with here. 

The appropriate ensemble of matrices is the circular orthogonal ensemble, but we will equivalently use the Gaussian 
orthogonal ensemble (GOE). Thus in the spirit of RMT (and statistical mechanics) we will equate any quantity with 
its ensemble average. Thus when Tr (Pa) = I = fN we get 

N-l 2-1 N-l 

< C >=<< Nf(l-f) T,Y,Y, l<n|<M| 2 |(n'|V> m >| 2 » • (5.1) 

* \ •* ' m=0 n=0 n'=l 

The double bracket represent average over the ensemble. GOE for eigenfunctions is particularly simple with the 
assumption that each vector covers an TV dimensional sphere uniformly. More details that are crucial to the following 
can be found for instance in Jl3| ] . 

Thus under the assumption that there are no favoured eigenfunctions we get 

< c >= « |<r# ro >| 2 |<n'|Vm>| 2 » (5.2) 
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We change to a different notation purely for the sake of convenience. We denote the the n 
vector A by x n \. Then using the result ]13| 



2 2 



» = 



1 



we get 



< c> 



N 



N 



N(N + 2) 



1 - 2/N. 



th component of the 



(5.3) 



(5.4) 



The final approximation is for large values of N (the classical limit). We see that this is consistent with our simulations 
for the average value of the fluctuations. It predicts that for generic systems the average is less than unity and that 
the leading correction is proportional to h. This agrees with our result for the standard map but the bakers map has 
a different behaviour as the leading correction was found to be of order This is consistent with what has been 

observed by others so far: that the quantum bakers map is somewhat "less random" than other models of quantum 
chaos @. 

The variance is calculated in a similar manner, but leads to higher order correlations that appear to be a pure 
"two- vector" result in comparison to that used for the average. Any two eigenvectors of the ensemble are correlated 
due to the requirement of orthonormality and the full derivation of the reduced densities is not known to the author 
at present. Nevertheless much can be derived about the variance. Concentrating on S (which is trivially related to 
the variance) we get 

4 



5=«2 



m2>mi n— no 



no+l — 1 



n\ip m2 )(ip mi \n) 



» 



(5.5) 



S = N(N-1) « 



5>„ 



» 



(5.6) 



We have made use of the fact that GOE eigenvectors are real and we also require in the above that A ^ (i. We have 
also dropped the no term (or set it equal to unity) as the statistical theory does not distinguish between eigenvectors. 
It is clear therefore that S is a measure that may be called "incomplete orthonormality" , while the average is related 
to "incomplete orthogonality" . 

The sum over / terms is equal (due to orthonormality) to a sum over N — I terms and therefore statistically we 
derive that 



S(l) = S(N-l) 



(5.7) 



where we have explicitly written the dependence of S on I. We must note that this is not an exact condition of 
symmetry (from detailed dynamics) as is apparent from the numerical results, but is exact within the framework of 
RMT. The deviations from this symmetry could therefore provide an interesting measure of deviation from statistical 
behaviour. 

Further expansion of the fourth power in the equation of S and book-keeping leads to the following fourth order 
polynomial for S. 



S(l) 



N(N — 1) 



= a t l + a 2 l(l - 1) + 6a 3 l(l - I) (I - 2) + a 4 l{l - !)(/ - 2) (I - 3) 



(5.8) 



and 



oi = « xf x x1j X », 
a 2 = 3 « x 2 iX x 2 jX x^x 2 jti » 

2 2 

< C < C X ^yX jyX i^X j iyX ^> , 

Q-4 ^"C. X^\X yXXip^X jy^Xii/X XifxXjoc ^>^> . 



4 << X iX Xj X Xi^Xj fl >>, 



(5.9) 
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These are truly two-vector averages ( i and j indices can be also thought of as denoting vectors rather than 
components). None of the indices must be the same in each of the above expressions. We see quite easily that 



a 4 = 3a 3 /{3- N). (5.10) 

This follows on summing over the index A in the expression for 04 and using orthonormality. Further we have that 

(N-l) « < v r>r, ; ,.r ;/ , »= - « xf x x% » . (5.11) 

This follows from summing over the index //. 

Further using the symmetry in the form of S gives 

ai + (N- l)o 2 + 3(iV - 1)(N - 2)a 3 = (5.12) 

Thus putting together this information leads to 



S(f = I IN) 
N(N-l) 



= /(!-/) 



2 



-N 2 (N- 1)04 



+ / 2 (1- f) 2 a 4 N 4 . (5.13) 



The symmetry about / = 1/2 is now obvious. a\ can be evaluated exactly using the methods in fll3|] and gives 

ai = 9/{N(N + 2)(N + 4)(iV + 6)) - 9/iV 4 . (5.14) 

The correlation 04 is unfortunately beyond the author, but it is sufficient to observe that 

a 4 — 3 << xaXjxXi^Xjn » 2 = 3/N e (5.15) 

The corrections to this are of order iV~ 7 . The factor 3 is neccessary as the number of different ways of choosing 
2 indices from 4 is 6 and choosing any two results in the automatic choice of the other pair. This new two-vector 
correlation is easy to evaluate ( |p^| ) and we get the above result. 
Collecting all the correlations we finally arrive at 

S(f = l/N) = 3/ 2 (l - f) 2 + OiN- 1 ) (5.16) 

Thus we have proved that S is only a function of / in the large N limit. This functional form should match those for 
time-reversal symmetric systems with no special symmetries. This is illustrated in Fig. 6, and the fit is satisfactory. 
The standard deviation of the fluctuations is for such systems thus given simply by 

a = V3/N + 0(N- 2 ), (5.17) 

and implies the independence of the width of the fluctuations from the size of the partitions. In contrast the numerical 
results show a strong / dependence for systems with special symmetries such as parity |s) . Further work is underway 
for elucidating features of the time-reversal breaking case as well. 



VI. SUMMARY 

We have studied correlations between operators in quantum chaotic systems. The relaxation towards an oscillatory 
state of decorrelation or "equilibrium" has been analysed numerically and also from the point of view of RMT. Certain 
universal features have been found and in the case of time-reversal symmetric systems these have been derived using 
RMT. 

There are many natural extensions to this work, the most important being consideration of general density operators 
rather than projectors. Further it has been found in B that the relaxation process is Gaussian implying that the 
mean and standard deviation characterise the process completely. The question of energy has not arisen in this work 







because of our models. Therefore the generalization of this work to Hamiltonian flows with infinite Hilbert spaces 
is of considerable importance. Connections have to be made between these approaches in quantum chaos and the 
existing methods of quantum statistical mechanics. 

The effect of scarring of wavefunctions by classical periodic orbits tends to localize many quantum states. It would 
be of interest to see if the measures defined herein are sensitive to this phenomenon. We have already noted deviations 
from RMT in the case of the quantized bakers map which is known to have a large number of scarred eigenstates. 
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FIG. 1. The relaxation fluctuations as a function of time for the symmetric standard map with K — 20 and (a) N = 100 (b) 
N = 200. 
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FIG. 2. The average of the fluctuations for (a) the quantum bakers map quantizing the (2/3,1/3) Bernoulli scheme and (b) 
the standard map with K — 20. Shown is the logarithm of the deviation from unity as a function of log(iV). 
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FIG. 3. The scaled standard deviation aN as a function of the kick strength for (a) the symmetric standard map, (b) the 
symmetric Harper map. The solid line corresponds to N = 300, the dotted line to N=200 and the dashed line to N=100. 
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FIG. 5. The eigenfunction sum S for the Harper map with 3 = 8. The symmetric Harper map with (o = 1/2, b = 1/2) is the 
solid line while the case (a — 0.35, b = 1/2) is the time-reversal symmetric case with no special phase space symmetries and 
the dotted line (a = 0.35, b = 0.35) breaks also time-reversal symmetry. In all cases N — 200. 
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FIG. 6. Comparison of the prediction from RMT, the solid line, with the asymmetric standard map 
(N=200,K=20,a=0.35,b=0), the crosses. 



13 



